Relativistic corrections to the Pionium Lifetime 
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Abstract: Next to leading order contributions to the pionium lifetime are considered within 
non-relativistic effective field theory. A more precise determination of the coupling constants is 
then needed in order to be consistent with the relativistic 7r — -k scattering amplitude which can 
be obtained from chiral perturbation theory. The relativistic correction is found to be 4.1% and 
corresponds simply to a more accurate value for the non-relativistic decay momentum. 

PACS numbers: 03.65.N, 11.10.S, 12.39.Fe 

In the DIRAC experiment which is underway at CERN, one plans to measure the 
pionium lifetime with an accuracy of 10% or better |T[[Q]. The dominant decay proceeds 
through the strong annihilation 7r + + 7r~ — ► 7r° + 7r°. Since the momenta of the final 
state particles is given by the mass difference Am = m + — mo between the charged and 
neutral pions, the process is strongly non-relativistic. In lowest order the decay rate follows 
directly from the corresponding non-relativistic scattering amplitude which can be written 
as 



in the center-of-mass frame where the energies of the pions are Eq = E + = m + + p 2 /2m + 
when p is the momentum of the charged pions. The S-wave scattering length a and 
the slope parameter b include both higher order chiral corrections and isospin-violating 
effects from the quark mass difference m u — m^ and short-range electromagnetic effects. 
At threshold the momentum p = and the full scattering amplitude with long-range 
Coulomb interactions removed is then given by just this scattering length. It can be 
written in terms of the more conventional isospin-symmetric scattering lengths ao and a2 
in the isospin 1 = and isospin 1 = 2 channels as a = ao — 02 + Aa where Aa includes 
these symmetry-breaking effects. 

Since the charged pions in pionium are supposed to be bound in a 15 Coulomb state 
\P(r) with relative momentum 7 = am + /2, the annihilation takes place essentially at 
threshold. In lowest order the transition rate is then given by just the scattering length a, 
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where |\&(0)| 2 = j 3 /ir gives the probability to find the two particles at the same point. A 
measurement of the pionium lifetime will thus give an experimental determination of the 
scattering length. For this to be meaningful, corrections to this lowest order formula must 
be calculated so that the theoretical lifetime has an uncertainty substantially smaller than 
the one in the experiment. 

In a recent paper we have described this system in the new framework of non- 
relativistic effective field theory Q. This approach offers a compact and more systematic 
approach to obtaining the many different higher order corrections to the pionium life- 
time compared with previous methods which to a large extent were based on covariant 
methods 0. As shown by Holstein||] to lowest order in the interactions this approach is 
equivalent to just using effective couplings in non-relativistic quantum mechanics. 

The effective field theory for non-relativistic pions is based upon the free Schrodinger 
Lagrangian Cq = ir<(idt + V 2 /2m 7r )-7r which corresponds to the propagator 

for a particle with energy E and momentum k. The interactions we consider are contained 
in the Lagrangian) 2 ] 

Li n t = -^Co(vr+vr* 7r 7r ) + ^C 2 {-n* + ^-noV 2 -kq + ir* + V 2 tt1tt tt ) + h.c. (4) 

where the gradient is defined as V = l/2(v — v). It makes it possible to calculate the 
scattering amplitude for ir + + 7r~ — » ir° + 7r° to order p 2 which then must agree with the 
definition (jl|) to this order. The result of this matching for the first coupling constant is 
then 
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when we neglect smaller rescattering corrections Q. Similarly, one finds for the derivative 
coupling C2 = (87r/3m^_)(6 — a). These values are more accurate than the ones used 
previously^]. 

With the value of the coupling Co now determined, one can calculate the second order 
correction AE to the ground state energy of pionium from the bound state diagram in 
Fig.l as first pointed out by Labelle and Buckley ||. It is found to be imaginary with a 
resulting decay rate of T = — 2ImAi?. To lowest order in the mass difference Am and 
ignoring the small binding energy, this gives the zero-order result (0). 

Since the pion scattering lengths are set by the natural size the counting rules 

needed to estimate the magnitudes of contributions appearing in different orders of pertur- 
bations theory, are simple. The energy E in the propagator will be of the order Q 2 when 

2 We define here Co with opposite sign to what we used in H 
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Figure 1: Next- to- leading order contribution to the ground state energy level shift which gives 
the leading order decay rate. The filled circle at the vertices denotes the Co coupling. 



the characteristic momentum in the process is Q. As a result, the propagator scales as 
1/Q 2 . For the same reason the four-dimensional volume integration J d 4 k will scale as Q 5 . 
The loop diagram in Fig. 1 thus scales as Q since it involves two one-particle propagators. 
This is the leading order contribution to the decay rate. 

To next order in the effective theory we must include the contribution from the deriva- 
tive coupling C2 in (||). It will follow from the diagram in Fig. 2 and is seen to scale as Q 3 
again using dimensional regularization^]. The contributions to the decay rate from these 
two diagrams are thus found to be 

r = ^|*(0)|V2Amm (C 2 + 2C C 2 Ammo) (6) 

With the above value for the two coupling constants we then simply get the lowest order 
rate @. The correction due to the slope parameter || thus disappears with the more 
precise values for the matched coupling constants used here. It results from a cancellation 
between a next-to-leading order p 2 contribution and a Am contribution which in this 
particular process are of the same order. 




Figure 2: Next-to- next-leading order correction to the ground state energy level shift which gives 
the next-to-leading order decay rate. The filled square at the left vertex denotes the Ci coupling. 



To this order in perturbation theory relativistic effects must also be included. These 
have previously also been considered in more covariant approaches 0]. Here they will arise 
from the lowest order relativistic correction to the free Lagrangian Co which now should 
be taken to be 

J V 2 V 4 \ 

This new interaction will modify the propagators in the bubble of Fig.l as shown in Fig. 3. 
Since the diagram now involves three propagators, it scales as Q 5 Q 4 Q~ G = Q 3 . This is 
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Figure 3: Relativistic correction to the decay rate. The cross denotes the relativistic coupling 



also confirmed by an actual evaluation of the diagram which is given by 

Here ^(p) is the Fourier transform of the bound state wavefunction. The integrations 
over momenta p and q now give just ^(O)) 2 . Using dimensional regularization for the 
divergent integral over the loop momentum k, we obtain the finite result 

AE (rel) = _^0 Am ^ Ammo |^( )|2 (g) 
327T 



In terms of the decay rate, it corresponds to 

Ar (reO ^ 5Am 
T 4m+ 
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and amounts to 4.1%. By its very nature it can also be derived from using more covariant 
methods^. Finally, the same relativistic interaction will also act on the external legs of 
the diagram in Fig.l. But the resulting correction is then of the order a 2 and can thus be 
neglected here. 

The relativistic correction (|l(]) follows also directly from the available phase space for 
the 7r°7r° final state since there is no energy dependence in the annihilation amplitude to 
lowest order. Each 7r° has the energy Eq = mo + k 2 /2mo — k 4 /8m,Q when we include the 
next-to-leading order term in the momentum expansion. The decay rate will then involve 
the integral 

f d 3 k 

r oc J ^ 5{2m + - 2E ) (11) 

when we ignore the small binding energy. The argument of the delta-function will now 
have two zeros of which one represents an unphysical, high-momentum state. Keeping 
only the contribution from the physical state, we recover exactly the additional term (JToj) . 

In a recent paper by Gall, Gasser, Lyubovitskij and Rusetsky0 the pionium lifetime 
is also calculated from the non-relativistic Lagrangian used above. But instead of using 
bound state perturbation theory based upon the standard Coulomb wavefunction as done 
here, they determine the properties of the bound state by calculating the complex pole 
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on the second Riemann sheet of the corresponding T-matrix. In this way they derive 
higher order corrections which are not considered here. Our results are consistent with 
their general form of the decay rate where the relativistic correction (|To| ) is seen to be the 
first term in the expansion of their decay momentum p*. The higher order corrections 
they have derived should follow in the present approach from including the Coulomb- 
interactions between the charged pions and rescattering effects corresponding to extra 
bubbles in the diagram Fig.l. 

Almost at the same time another calculation of the lifetime was completed by Eiras 
and Soto[||. This is again based upon the same effective Lagrangian which is now further 
reduced by integrating out degrees of freedom with momenta of the order \[2Kmm§. 
It then becomes an effective theory for only charged pions with contact and Coulomb 
interactions. Their result for the lifetime is also in agreement with what we have obtained 
here. 
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